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ABSTRACT 



An investigation of the dynamic stability of a typical fixed wing 
VTOL aircraft operating in the transition region of flight was conducted 
by means of analog simulation. The Ryan Aeronautical Company XV-5A 
Lift Fan Research Aircraft was used as a basis for the analysis. Wind 
tunnel tests on a simple 1:24 scale model of the aircraft provided static 
aerodynamic data throughout the angle-of-attack range -90° < o< w < 
+90°. Non-dimensional stability derivative data was extended through- 
out the angle-of-attack range. Analog simulation at relative velocities 
of 100, 60, and 20 knots was conducted using dimensional parameters. 
Recordings of the simulation runs provided period, frequency, damping, 
and dutch roll parameter information which was analyzed to determine 
areas of flight instability. The variations of these parameters with both 
angle-of-attack and velocity is discussed. The need for automatic 
stabilization is indicated. 
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1.0 INTRODUCTION 



The versatility of a fixed wing VTOL aircraft will be enhanced 
if it has the capability of flying at arbitrary combinations of vertical 
and horizontal speeds in the transition from hover to conventional for- 
ward flight and vice versa. This implies flight at arbitrary angle-of- 
attack in the range - 90°^ o< 90° and at arbitrary relative air- 

speeds. Throughout this flight envelope wide variations of aerodynamic 
characteristics can be anticipated. The dynamic behavior of the air- 
craft, therefore, can be expected to change radically. 

The purpose of this analysis is to investigate the dynamic 
stability of a typical fixed wing VTOL aircraft in this transition region 
of flight throughout both the angle -of-attack range - 90°.5c>< w <+90 o and 
the speed range 20 < 100 knots. The dynamic response of such an 

aircraft after being perturbed from a reference flight condition is 
analyzed in an effort to establish the effects of the varying aerodynamic 
forces and moments on the behavior of the aircraft. Areas of insta- 
bility as well as the degree of stability of instability as functions of 
both speed and angle -of-attack are obtained. The feasibility of manual 
control in the high transition speed - high angle -of-attack region of flight 
is considered and the necessity for automatic control is indicated. 

The aircraft chosen as a basis for the analysis is the Ryan 
Aeronautical Company XV-5A. This aircraft, built as a research 
aircraft for the U.S. Army, employs the General Electric X- 353-5 
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Life Fan Propulsion system. The aircraft is capable of both conven- 
tional flight at high subsonic speed and VTOL operations. It is not the 
intent of this analysis to define the dynamic behavior of the XV-5A per 
se but to indicate the dynamic behavior of a typical VTOL aircraft 
using the XV-5A as an example. For., this reason, although the XV- 
5A is equipped with manual controls, only stick fixed stability was 
considered. 

The analysis required extrapolation of stability derivative data 
throughout the angle -of-attack regime. To this end, static wind tunnel 
tests were performed to obtain the basic lift, drag, and pitching mo- 
ment data. The longitudinal and lateral-directional stability derivative 
data was then extended by means of theoretical analysis, empirical 
relations, and estimates to cover the specified angle -of-attack range. 
Finally, an analog simulation study, using dimensional parameters, 
was conducted over the ranges of angle -of-attack and flight velocities 
indicated previously. 

This investigation was conducted at the United States Naval 
Postgraduate School, Monterey, California, during the period August, 
1964, through April, 1965. 

The authors wish to express their appreciation for the advice 
and assistance given them by Professor E. J. Andrews of the Aero- 
nautical Engineering Curriculum. Appreciation is also extended to 
the Ryan Aeronautical Company for making available the XV-5A data 
used in the analysis. 
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2.0 DISCUSSION 



An investigation of the longitudinal and lateral-directional dynam- 
ic stability of the XV-5A VTOL aircraft in the transition flight regime 
was conducted by means of an analog simulation of the appropriate equa- 
tions of motion. The aircraft was taken to be in the transition configu- 
ration with 45° flap deflection and 15° aileron droop. Longitudinal po- 
sition of the center of gravity was taken to be at Station 246. 0. The 
flight condition was assumed to be at sea level and out of ground effect. 

In the performance of a transition maneuver from hover to con- 
ventional flight, and vice versa, a fixed wing VTOL aircraft may en- 
counter the complete range of angle -of -attack -90°^o< w ^ +90° due to 
the infinite number of possible combinations of vertical and forward 
speeds which might be achieved. In the present analog study of the dy- 
namic behavior of the XV-5A, the responses to a step or impulse per- 
turbation of pitch attitude and lateral velocity were measured for spe- 
cific steady state conditions of angle -of -attack and relative velocity. 

Time permitted a detailed investigation of only the positive angle-of- 
attack range +90° corresponding to a descending transition 

from conventional flight to the hover. Some qualitative conclusions 
were drawn, however, concerning the dynamic behavior in the negative 
angle -of -attack range. 

Only the stick fixed dynamic stability of the XV-5A aircraft was 

investigated. Further, it was assumed that the symmetric and assym- 
metric motions of the aircraft were uncoupled, allowing separate 
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analyses of the longitudinal and lateral-directional dynamics. 

Before the analog simulation could be conducted, it was neces- 
sary to obtain estimates of the aerodynamic characteristics of the air- 
craft over the complete angle -of -attack range -90°^ <x w s +90°. Reli- 
able aerodynamic data was available for the low angle of attack range, 
herein considered to extend from 0^=0° to ©<^=16 0 . (12, 13) For the 
range of large angles -of -attack, defined to extend from C »< V =16 0 to ©< w 
=90°, the available low angle -of -attack data was extended by various 
methods described in detail in following sections. Estimates of some 
of the longitudinal stability derivatives were obtained from lift, drag, 
and pitching moment data resulting from static wind tunnel testing of a 
very simple model of the XV-5A aircraft. The wind tunnel tests of this 
model at large angles -of attack are discussed in Section 2.1. Estimates 
of the non-dimensional stability derivatives, presented in Section 2.2, 
were based on a system of stability axes, or wind axes, since the avail- 
able data for low angles -of -attack was expressed in this system. It was 
desirable, however, to use dimensional stability derivates referred to 
body axes in the analog simulation. The stick fixed, uncoupled equa- 
tions of motion, referred to body axes, are defined in Section 5.1. The 
relationships for converting the various non-dimensional derivatives to 
dimensional form and the equations for transforming these derivatives 
from stability axes to body axes are given in Section 5. 3. Also includ- 
ed in this section is a digital computer FORTRAN program designed to 
perform these conversions. The analog circuits for the simulation of 
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the equations of motion are described in Section 5.4. The analysis of 
the results obtained is presented in Section 2. 3. 

Although the phugoid oscillation is usually assumed to occur at 
constant angle -of -attack in conventional stability analyses, it was not 
clear at the outset of this investigation that this would be the case for 
flight at angles -of -attack above the stall. Since most of the longitudi- 
nal stability derivatives vary with in a highly nonlinear manner 

*r 

over some portions of the angle-of-attack range, it was considered 
necessary to simulate the variation of these coefficients by the use of 
variable base diode function generators in the analog Computer circuit. 
The variations of the longitudinal dimensional derivatives for a devia- 
tion of angle-of-attack of jf 1° from a given reference angle were pro- 
grammed into the function generators. 

2. 1 Wind Tunnel Tests 

In order to obtain the basic lift, drag, and moment data of the 
XV-5A aircraft throughout the angle-of-attack range -90°:<cx^< +90°, 
static wind tunnel tests were conducted on a scale model of the aircraft 
The model, shown in Fig. 2. 1, was constructed to a scale of one inch 
equal to two feet using data obtained from RAC. (12) A three view 
drawing of the XV-5A aircraft is shown in Fig. 2.2. No attempt was 
made to construct a sophisticated reproduction of the aircraft as far as 
airfoil sections, flaps, and controls were concerned. Rather the em- 
phasis was placed on obtaining the proper planform, scaled dimensions 
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and streamlined contours. It was anticipated that above the stall angle- 
of -attack the lift, drag, and pitching moment produced on the aircraft 
would be the result of a 'flat plate' force effect since flow over all air- 
foil components would be completely separated. Planform would there- 
fore be most important. Below the stall angle -of -attack it was antici- 
pated that the data of RAC would be used in the analysis. (12) Never- 
theless, data was obtained in the region of low angles -of -attack in order 
to compare results with RAC as an indication of the adequacy of the 
model used. 

Model testing was conducted in the Aerolab Subsonic Wind Tun- 
nel located in Building 225 at the U. S. Naval Postgraduate School. A 
planform layout of the tunnel is shown in Fig. 2.3. A view of the tun- 
nel test set-up is given in Fig. 2.4 showing the complete model mounted 
by a three point suspension at a positive angle -of attack of approximate- 
ly 45°. Various sting mounts other than the one shown were used in the 
lower and negative angles -of-attack ranges. The mechanical balance 
used to measure the forces and moments is also shown. 

Model test data was obtained in both the WBVH and WBV con- 
figurations in the angle of attack ranges of 0° <cx w ■< +82° and -82° "S. 

S -43°. Due to the physical constraints of the mounting system 
employed, angles of attack above 82° in magnitude were unattainable. 
Data in the range -43°<o< w < 0° was not obtained because of sched- 
uling difficulties of the tunnel facilities. 

Initial testing was conducted in the angle -of-attack range 
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42°<o< w < 82° at three flow velocities: 30, 60, and 120 feet per 

second. Preliminary analysis of the data obtained showed no notice- 
able variation of the results between the various flow velocities. The 
remaining tests were conducted, therefore, at 120 FPS. 

All tests were conducted in a clean aircraft configuration. 
Results of the tests are presented in non-dimensional form in 
Figs. 2. 5 and 2. 6. Fig. 2. 5 shows Cj^ and C-^ as functions of wing 

angle -of -attack, o< w , while Fig. 2.6 shows C m as a function of o< . 

© 

Both figures include RAC data for comparison. (12) 

In view of the relative crudeness of the model used in the tests, 
agreement with RAC data is considered good, particularly for C m 

c § 

The lift curve slope is in very good agreement but wing stall occurs at 
lower values of o< and C T . The C n data is most in error since 
no attempt was made to correct for flow into the engine intake duct. 
Consequently, the duct of the model presented a flat plate surface to 
the flow which produced a rather severe drag increase in the low angle- 
of-attack region. 

Above the stall angle -of -attack, both C^ and C^ drop off 
abruptly indicating flow separation and a breakdown of the circulation 
lift mechanism. This decrease continues to an angle-of -attack of ap- 
proximately 2 0° when both C-^ and Cj^ begin to rise once again. How- 
ever, above the stall angle -of -attack lift and drag forces are predom- 
inately force components, along specified axes, of the total force on 
the aircraft resulting from a 'flat plate drag' effect. This is borne 
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out by the fact that Cl and Cq are equal at approximately c>< ^ =+46° 
and Cl reaches a maximum value, above the stall, at approximately 
o^=+47°. A pure ’flat plate' effect would produce these effects at 
o<' VsJ =+45°. To be noted also is the high value of Cl reached in this region. 
The value of Cl equal to 0.9 at c*' VJ =+47 0 is only approximately ten per- 
cent lower than CL max ^h e stall where the wing is producing lift 'ef- 
ficiently 1 . As o<^approache s +90°, Cl approaches zero indicating 
the flat plate nature of the results. 

In the angle -of-attack range + 36°<o<^ <+40° CLWBV * s c °i n_ 
cident with CLwBVH indicating a complete blanking of the horizontal 
tail. On either side of this angle-of-attack range, the ClwbVH an< ^ 
CLWBV curves diverge. As o' approaches +90°, however, CLwBV 
approaches zero as did CLWBVH • The appearance of the deep stall 
region where the horizontal tail is completely ineffective is equally ob- 
vious in both the Cq and C mc g plots. In this region of angle-of-attack 
it may be expected that both the static and dynamic stability of the air- 
craft would deteriorate sharply. The static stability deterioration is 
most vividly shown in the C m cg plot. 

Although data in the angle-of-attack region - 43°^o< w 0° 

was not obtained, it can be expected that the deep stall region would 
not appear in this range since the horizontal tail would not be blanked 
at any negative angle-of-attack. This result is indicated by the avail- 
able data in that ^LwBVH 9 ^DWBVH > an< ^ ^ m cg WBVH are s i - 
cantly different from CLWBV’ ^DWBV» and Cm C g WBV in the re gion 
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of <=*< = - 43° and show no tendency toward coincidence. Once again 
Cl and Cq become equal at cx^ = - 44° with Cl indicating a peaking 
tendency at approximately °< w = - 43°. 

Finally, above the stall angle -of-attack, it can be seen that the 
slope of the C mc g curve is such that a high degree of static stability is 
inherent in the aircraft throughout both the positive and negative angles - 
of-attack except in the deep stall region of positive angles -of-attack as 
previously discussed. 

With this data available it was possible to predict the static 
aerodynamic characteristics of the aircraft in the transition configura- 
tion throughout the angle -of-attack range -90°^ 5 +90°. The final 

CLWBVH an d CDWBVH curves used in this analysis are presented in 
Fig. 2. 7 while the final Cm C g yfgyjq curve used is shown in Fig. 2. 8. 
These curves are made up, in general, of three distinct parts as indi- 
cated. For angles-of-attack +17°, RAC data is used. (12) 

For angles-of-attack + 1 7°< < +82° and -83°^<»<' w ~ -43° the wind 

tunnel data just discussed is used. For the remaining angle -of-attack 
regions, estimated variations have been used. These areas are pri- 
marily the negative angles-of-attack from 0° to -43° and that portion 
of the lift curve immediately following the positive stall angle -of-attack. 
Due to the 'flat plate 1 effect above the stall, the presence of flaps and 
drooped ailerons will have only a minor effect on the Cl an d Cq data 
obtained from the wind tunnel tests without these features incorporated 
on the model. For this reason, the data from these tests was used 
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directly for the aircraft in the transition configuration. 

Fig. 2.9 presents the variation of CLt and Cp)t> based on hori- 
zontal tail area, as a function of angle -of-attack. The variations for 
angles -of-attack +17°^ o< < +82° and -83° ^ «=< <* -43° were ob- 

tained from the differences (Cl wb y H - CLWB v) and( CDWB VH“ c DwB 
of Fig. 2. 5 and converted to a basis of tail area. These portions of the 
plots are presented as a function of wing angle -of-attack, , while 
that portion below the stall, obtained from RAC data, is presented as a 
function of the tail angle -of-attack, . (12) Use of the data in the latter 

region for prediction of the effect on the aircraft as a whole requires the 
standard corrections for and to establish the actual angle- 

of-attack and Cl of the tail. Data above the stall angle -of-attack may be 
applied directly to the aircraft as a whole since tail blanking effects are 
obviously inherently included in the presentation. 

As in the case of the complete model, CLt and Cp)t are equal at 
angles -of-attack of approximately +46° and - 43°. As <=< approaches 
90° CLt approaches zero while Cp|- approaches a peak value. Above the 
stall angle-of-attack, CLt reaches a peak value at angles -of-attack of 
-43° and +70°. The tail blanking effect in the positive angle-of-attack 
range results in a lower peak value of C-^. than in the negative angle-of- 
attack range and at an angle far removed from the expected 1 flat plate 1 
45°. 

The composite curves of Fig. 2.9 are those used in the 
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following parts of the analysis. 



2.2 Non-Dimensional Stability Derivatives 

As a prerequisite to the dynamic stability analysis, it was 
necessary to estimate the magnitude of each of the relevant non- 
dimensional stability derivatives over the range of angle -of-attack 
from 0° to 90°. Good estimates of all of the non-dimensional sta- 
bility derivatives were available for positive low angles -of-attack 
up to the stall ato^=l6° (12, 13). Various methods, described in 
this section, were used to extend this data for large angle -of-attack. 

It should be noted that the lateral-directional rate derivatives are 
the most difficult to predict, either analytically or experimentally, 
even at low angles -of-attack. The complex empirical techniques 
used to estimate these derivatives below the stall are valid only in the 
range of angle -of-attack for which the slope of the lift curve, C L » 
is constant (2). In the absence of accurate means of predicting 
the values of the lateral directional rate derivatives at large angles - 
of attack, it was necessary to use elementary methods of analysis to* 
indicate their approximate variation over this range. 

2.2. 1 Longitudinal Non-Dimensional Stability Derivatives 

2.2. 1. 1 Static Longitudinal Derivatives, c, , c , c 

°c* ^ 

Values of the derivatives C , C , and C for the WBHV and 

i-t* o«*’ 

WBV configurations were obtained from the lift, drag, and moment 
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curves described in the previous section. These derivatives are 



plotted in Figs. 2. 10 through 2. 12 inclusive for the positive angle- 
of-attack range, 0° £ ^90°. Also included in Fig. 2. 10 is 



2.2. 1.2 0 } tt\ © 

L <* <=< 

The lift and pitching moment coefficients due to rate of change 
of angle-of-attack arise from the lag in downwash at the horizontal 

tail. Equations (2- 1) and (2-2), used in estimating C L o and o are 

°< °< 

derived in Section 5. 2. 



c L . =ec L 'w |± 4 |^- 

L e* L o<-t <-t S w C d°< 






-c h. 

c 



(2-1) 

( 2 - 2 ) 






It has been shown that the slight variation of with angle -of-attack 

is approximately linear below the stall for the XV-5A. (12) Since down- 
wash results from circulatory flow about the wing and since the wind 
tunnel tests, discussed in the previous section, showed that circulation 
decreased rapidly above the stall, it was concluded that C. 0 and 



'ex' 



approached zero as angle-of-attack increased above the stall. The 



data of Ryan for C and C, were therefore extended to reach zero 
at ^28° and were replotted in Figs. 2.13 and 2.14. (13) Both of 



these derivatives were assumed to be qualitatively similar at positive 
and negative angles -of-attack. 
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2. 2. 1 . 3 C > C ^jq 

The wing and horizontal tail make significant contributions 
to the derivatives C, and C . Curvature of the flow relative to a 
translating wing occurs when a pitch rate, , is imposed on the 
wing. This effect can be thought of as an equivalent change in cam- 
ber which alters the lift and pitching moment characteristics of the 
wing. The empirical relation of Equation (2-3) was used in estim- 
ating this effect (6). 

The quantity varies over the angle -of-attack range according 

to the geometric construction of Fig. 2. 15. 

The contribution of the horizontal tail to the derivative results 
from the change in angle -of-attack due to the vertical velocity com- 
ponent of the tail and is given by Equation (2-4) which is derived in 



Section 5. 2. 



>f ^ 



r - 2C 17 — — 

% s w c 



(2-4) 



The variation of with <=*<^is shown in Fig. 2. 16. The sum of the 

components and C were plotted in Fig. 2.17 over the positive 

+\n 

and negative range of angle -of-attack. For negative angles -of-attack 
the horizontal tail does not pass through a region of deep stall or 
separated flow from the wing. The dotted portion of the curve of 
Fig. 2. 17 represents an estimate of C^ Q in the range of negative 



angle -of-attack corresponding to the positive angle deep stall range. 
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The contribution of the wing to the pitching moment coefficient 



due to pitch rate, , was estimated by means of the empirical Equa- 
tion (2-5). (6) 






= -c L . 



W 



p ! 


[A+eCosA| 


1. 


Lb‘ 


r A j 


I 1 






The inboard and outboard sections of the wing of the XV- 5A have 

* 

different sweep angles. A mean sweep angle, -A = 20°, was assumed 
for use in the above equation. 

The effect of the horizontal tail on C was estimated from 
Equation (2-6). 



c - _ r M 
St £ 



(2-6) 



The derivative C is plotted in Fig. 2. 18 against angle-of-attack. 
As in the case of C , C is not altered by the deep stall condition 
in the negative angle-of-attack range, as indicated by the dotted 
portion of the curve. 



2.2. 1.4 C , C , 

The non-dimensional derivatives of forces normal and paral- 
lel to the wind axis and of the pitching moment due to airspeed are 
given by Equations (2-7), (2-8) and (2-9) which are derived in Section 

5.2. 



C 



-£C 



o 



(2-7) 
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( 2 - 8 ) 





(2-9) 



Since the aircraft is assumed in trimmed flight with thrust pertur- 



tion of Equations (2-7, (2-8) and (2-9) were Mach number effects 
which are of no interest in this study. 

2. 2. 2 Lateral-Directional Non-dimensional Stability Derivatives 



eral angle -of-attack, @ , on the fuselage and on the vertical tail. 

The wing was assumed to provide a negligible contribution to this 
derivative. At low angles -of-attack a large resultant negative side- 
force, acting close to the nose of the fuselage, can be expected when 
a sideslip angle is established. AsC>< is increased, the effective 

w 

sideslip decreases, the center of pressure of the sideforce will 
move aft on the body, and the magnitude of the sideforce will decrease. 
In the limit, as <=><^approaches 90°, only a small sideforce, acting at 
approximately the geometric centre of the body, will be exerted. For 



ly determined curves. (12) This data was extropolated to an estim- 



bations ignored, C was neglected. Also neglected in the deriva- 



2 . 2 . 2 . 1 




The sideforce due to sideslip results from imposing a lat- 



angles-of-attack below the stall, C 




obtained from exper imental- 



ated small value of C at °< = 90° as shown in Fig 2. 19. 

\ w 
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A negative sideforce is exerted on the vertical tail by a side 



slip. A slight variation of S, with °< , indicated in the experimen- 
tal results for low^, is probably caused by fuselage interference. 
(12) As angle-of-attack increases to the stall, part of the vertical 
tail becomes immersed in the separated flow behind the wing. The 
effective tail area and C are decreased. This deterioration be- 

comes more severe as o< is further increased and more of the tail 

w 

area is covered by the separated flow region. Since it appears that 
at least part of the vertical tail will be in unseparated flow at all 
times, it was estimated that C would decrease to a minimum of 

i|p v 

15 percent of its average value for low angle s -of-attack. It was 
noted that the vertical tail would have minimum effectiveness at a 
smaller angle-of-attack than would the horizontal tail. For this for- 
mer, this angle was estimated to be 32°. As the vertical tail em- 
erges from beneath the wing, at very high angles -of-attack, it tends 
to regain effectiveness. However, for very large wing angle-of- 
attack, the effective lateral angle-of-attack produced by a given side- 
slip velocity will be very much greater than the stall angle of the 

vertical tail. Therefore C u will be very small for large ; 

w 

Ou at = 45° was therefore estimated to be about 35 percent of 

decreased 

Pv 

toward zero at o< = 90°. 

W 

The estimated variation of CV witho< is plotted in Fig. 2. 19, 

w 

based on available data for C u at low angle-of-attack. (12) 



its magnitude below the stall. 



Above C>< = 45°, 0 o 

w il 
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2. 2. 2. 2 C» 



The yawing moment due to sideslip comprises contributions 
from the fuselage and vertical tail. These moments are the products 
of the sideforces described above acting on lever arms about the air- 
craft centre of gravity. The moment arm of the fuselage sideforce, 
measured from the fuselage center of pressure aft to the centre of 

gravity, decreases to zero asO< approaches 90°. The moment arm 

w 

of the vertical tail sideforce is measured from the aerodynamic cen- 
ter of the mean aerodynamic chord of the vertical tail. Due to the 
location of this point relative to the aircraft center of gravity, the 

vertical tail moment arm will not vanish at o< = 90°, but will be re- 
vs/ 

duced to a minimum of about 20 percent of the arm at C< = 0°. 

w 

The wing contribution to was considered to be negligible. 



Fig. 2. 20 shows the extrapolation of the variation of G at low 
angle -of-attack to the large angle range (12). 



2. 2. 2 . 3 Co 

y 

The dihedral effect of the wing and the sideforce acting on the 
vertical tail comprise the rolling moment due to sideslip. Experi- 
mental results show thatGp increased negatively with angle-of- 
attack up to the stall, apparently due to fuselage interference. (12) 
Above the stall, an abrupt reduction, in differential lift on the wings 
would be expected for a given sideslip. The reduction in lift results 
from the decrease in effective angle -of-attack change due to dihedral. 
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The component, Cg , is the product of C and the moment 
arm measured perpendicular to the relative wind axis, through the 
center of gravity, from this axis to the aerodynamic center of the 
vertical tail. For small angles -of-attack the aerodynamic center of 
the tail is above the relative wind axis and a negative rolling mo- 
ment results from a positive sideslip. For large positive*^ , 

w 

however, the position of the vertical tail is such that a positive 
rolling moment is caused by a positive sideslip. The length of the 
moment arm and therefore the value of Cn are zero at c< = 16°. 

Jr O W 



In the r 



ange 16°^ c>< <36°, C,p remains small due to blanking of 



the vertical tail by the wing. For very large 0< , although the side- 

W 

force on the vertical tail due to sideslip would be small, a signifi- 
cant rolling moment would be expected due to the very large mo- 
ment arm on which the sideforce acts. Note that the discussion 
refers to stability axes and that a rolling motion, about the X sta- 
bility axis at o< = 90° has become a yawing motion about the ^ 
body axis. 

For the negative range of©< , the moment arm will always be 

w 

such that Cp will be negative. No blanking of the vertical tail would 
be experienced in the range - 1 6°S - 36° as in the positive range. 

Rather, it would be expected that at very large negative angles -of- 
attack, the vertical tail w^.uld become pxfogre s sively blanked by the 
horizontal tail. 

The available data for C at low angle -of-attack was extended 
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to large angles by means of the above methods and the resulting 
curve is presented in Fig. 2.21 (12). 



2. 2. 2. 4 



C 



Mr 



The major contribution to the sideforce due to yaw rate,C u , 
comes from the vertical tail. The empirical methods mentioned 
previously were used to estimate the very small contribution of the 
wing to the derivative C^at low angles -of-attack. (13) The effects 

of wing and fuselage on were assumed to be negligible at large 

* )r 

angles -of-attack. For the range -90°^o< <90°, Equation (2-10) 

w 

was used to estimate Cy . (7) 

HI 



Mr v 



c - ec 



(2-10) 



Fig. 2. 22 shows the variation of 0 U with oC* w 

W 



2 . 2 . 2 . 5 C 






The yaw damping derivative, results from differential 

drag due to differential velocity over the wings and from the yawing 
moment due to the lateral angle -of-attack change associated with the 
yaw rate on the vertical tail. Equations (2-11) and (2-12), derived 
in Section 5.2 for assumed constant span wise lift and drag distri- 
bution, were used to estimate the wing and tail contributions to 
at large angle s -of-attack 



0 ^ = ~ Cp (0.8) 



(2-11) 



w 
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( 2 - 12 ) 



c 



n 






£ C 
c 



The variation of 




with c< w is 



shown in Fig. 



2. 23. 



2. 2. 2. 6 Cjp r 



The rolling moment due to yaw rate generated by the wing 
results from differential lift due to velocity changes over the wings. 
Handbook methods were used by the manufacturer to account for 
geometric effects and non-constant span wise lift distribution in estim- 
ating (3p for small angle s -of-attack. (2, 13) For large o' , 
Equation (2-13), an approximate relationship, was used. 



Cje 



c L 






( 0 . 8 ) 



( 2 - 13 ) 



The derivation of Equation (2-13) is given in Section 5.2. 

Because of lack of circulation at large angles -of-attack, the 
assumption of constant lift distribution used in deriving Equation (2-13) 
was considered to be more accurate than below the stall where circula- 
tory effects exist. 

The vertical tail contribution to < -JL was computed from 
Equation (2-14). (7) 



- e 



^2£v Q 

C iJ(3 V 



Pv 



(2-14) 



Cj?, is plotted against for positive and negative angles in 

Fig. 2.24. 
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2. 2. 2. 7 C 






The only significant component of the sideforce derivative 
due to rolling moment comes from the vertical tail. Equation (2-15) 
was used to determine c *. • (?) 



= 2 + C 



^ c 



(2-15) 



Pv ' ~*Pv dig' 



The second term of Equation (2-15) describes the effect on the verti- 
cal tail of sidewash due to roll rate. As angle-of-attack was in- 

a<r 

creased above the stall, A-pb was assumed to approach zero rapid- 

S fUo 

ly. At lowO<^there exists a small sideforce due to roll rate which 
arises from wing dihedral. This effect is very small compared to 

the vertical tail contribution and since it decreased with c< it was 

w 

neglected. C is plotted in Fig. 2.25 for the range - 90° ^ «=< ^ 

y-p w 

900. 



2 . 2. 2. 8 Chf 

Both the wing and the vertical tail have important effects on 
the yawing moment due to roll rate. The contribution of the wing 
comes from the change in magnitude and inclination of the resultant 
force vector on each wing. From elementary analysis considering 
constant span wise life and drag distribution, Equation (2-16) was 
obtained. 

r - (c - C ] ] ( 2 -i< 

c h fSy - )*, bro ) 

Derivation of this expression is shown in Section 5.4. For very low 
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angles -of-attack* this elementary analysis is in reasonable agree- 
ment with experimental data. (13)' At very high angles -of-attack, 
where circulation no longer exists, the above expression is also 

considered valid. It is important to note that C must be obtained 

°<x 

experimentally for high angles -of-attack since induced drag results 
from circulatory lift only. Therefore results of wind tunnel tests 
previously discussed were used in evaluation of C. at high angles- 
of-attack. The method just described predicts negative values of 
for angles -of-attack just above the stall where some circulation still 
exists. Experimental results for wings of similar planform indicate 
that C h will be positive in this region in contradiction of the theory. 

(8) Therefore from the data of this reference, C K . was estimated 

n t>v 

for 16 ° < ^ 32°. This estimate was faired into the very high 

angle-of-attack estimate and into the low angle data previously 
mentioned. 

The vertical tail contribution was determined from Equation 
(2-17). (7) 



c = £ + c h 

h fv C ^(3 V ™ 9 fuc 



(2-17) 



The variation of Cf^with c^^is shown in Fig. 2.26 for - 90 c ^°<^~ 90 c 
2. 2. 2. 9 Cj^ 

Differential lift on the wings due to differential change of 
angle-of-attack provides the largest component of the roll damping 
derivative. ^J^^ D was es timated from Equation (2-18) for large 
angles - of-attack. 
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(2-18) 



=“( C L„ +C d) w (|t|) 

The derivation of Equation (2-18) is given in Section 5.2. The contri- 
bution of the vertical tail, given by Equation (2-19), was found to be 
negligible. (7) 



Cjt. 



‘~Pv 

was plotted in Fig. 



= 2 



C il{3v 



c 
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3c- 

2Uo 



(2-19) 



2.27 for the range -90° ^ — 90°. 



2. 3 Analysis of Analog Computer Results 



The longitudinal response of the aircraft to perturbations in 
lateral velocity were measured for the discrete reference flight con- 
ditions of steady state angles -of-attack and relative airspeeds listed 
in Tables IX and XIII in Section 5.4. All computer runs were made 
for a steady state pitch attitude of 0 O = 10° except for one run at 

= 8°, V Q = 100 knpts, and © c = 0°. In the longitudinal analysis the 
perturbations U, , UT , o< , ©, and ^ were recorded on a multi-channel 
strip recorder. Perturbations of 'IT, , and Y were recorded 

in the lateral-directional tests. 

The first computer run of each of the longitudinal and lateral- 
directional series of tests was conducted for o< = 8°, V Q = 100 knots, 
and 0 O = 10°. These reference conditions correspond approximately 
to the minimum speed for which the period and damping parameters 
obtained from these initial computer runs were compared with the data 
of RAC for trimmed conventional flight at V Q = 94.6 knots and ©^—8. 1°. 
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For the longitudinal case, from Fig. 5. 12, the period of the 



phugoid oscillation was measured to be 39. 5 seconds, exactly the 
period given by the manufacturers report. (13) The measured time 
to damp to half amplitude was 20. 0 seconds for the initial run at 



first, resulted in a value of T 1/2 = 12. 3 seconds measured from 
Fig. 5. 13. These values of T \/2 are in reasonable agreement with 
the contractor value of T l/2 = 14. 6 seconds for 0 O = 8. 1°. 

Similarly, the initial lateral-directional test produced a 
period and time to damp to half amplitude for the dutch roll oscil- 
lation of 3.75 seconds and 2.5 seconds respectively, as measured 
from Fig. 5. 35, which compare very favorably with the RAC data 
of 3. 32 and 2.29 seconds respectively. (13) In addition, the dutch 



given value of .610. 

Once the analog circuits had been proven in the manner 
just described, the computer runs were carried out for the refer- 
ence flight conditions of interest. The results of each run are 
shown in Figs. 5. 12 through 5. 34 for the longitudinal mode and 
Figs. 5. 35 through 5. 57 for the lateral-directional mode. 



e. 



= 10°. A second run with 0 = 0°, but otherwise similar to the 

© 




was measured to be . 6 12 compared with the 



24 



2. 3. 1 Longitudinal Mode 



One of the most significant conclusions which can be 
drawn immediately from an examination of the longitudinal results 
is that the deviation of the angle -of-attack from its reference 
flight condition is very small for all runs. The variation is so 
slight that it is concluded that the longitudinal stability derivatives 
could be considered constant at any given reference angle -of-attack. 
Hence, in a future simulation of this nature, the function generators 
could be omitted from the analog circuit. The stability derivatives 
would be inserted as simple potentiometer settings thus simplifying 
the circuit. 

It was found that for angle s -of-attack at which the longi- 
tudinal response was oscillatory, the period varies considerably as 
a function of both angle -of-attack and airspeed. Table I lists the 
period and frequency of the longitudinal response at the various 
reference flight conditions. This data is not presented graphically. 
Such a presentation would be misleading since ranges of angle -of- 
attack exist where the motion is aperiodic. From Table I, however, 
it is clear that the period becomes quite short in the vicinity of the 
stall, increasing at very high and low angles -of-attack. It may also 
be seen in Table I that at low angles -of-attack the period decreases 
slightly with reduced airspeed while at large angles the opposite is 
true . 
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Wide variations in the damping characteristics were ob- 
served over the ranges of angle -of-attack and airspeed. For the 
purpose of graphical presentation of these variations, a damping 
parameter l/T was defined. The quantity, T, refers to the time to 
half amplitude, in the case of a stable oscillation, or time to double 
amplitude, in the case of an unstable oscillation or pure divergence. 

T, or T 2 was measured for each flight condition from Figs. 5. 12 
through 5. 34. The parameter l/T plotted as a function of angle-of- 
attack is presented in Fig. 2.28. The value of l/T = 0 in this figure 
represents the condition of neutral stability or undamped oscillation. 

At low angles-of-attack it is clear from Fig. 2.28 that, 
as airspeed decreases, the damping, and hence stability, deterio- 
rates. At = 8° the aircraft is stable at 100 knots, neutrally sta- 
w 

ble at 60 knots, and unstable at 20 knots. At the stall there is a 
sharp increase in the rate of divergence for all airspeeds. It is 
important to note that, in contrast to the low angle -of-attack region, 
the rate of divergence is greater at 100 knots than at 20 knots for the 

range of o' just above the stall. This phenomenon occurs because of 

w 

the change in polarity of major longitudinal derivatives such as 
and M^near the stall. It is not possible to be more specific about 
this point because of the different manner in which these derivatives 
vary with angle -of-attack. 

The second region of rapid deterioration of stability, shown 
by Fig. 2.28, is in the vicinity of 28°. This deterioration results 
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from a combination of several factors* At the stall the derivative 
^becomes positive and remains so up to an angle -of-attack of 
approximately 34°. The derivative M^remains negative until an 
angle -of-attack of approximately 26°. It then becomes positive and 

combines with causing the rapid divergence indicated at o< = 28°, 

+ w 

At an angle -of-attack of approximately 40° M^again becomes nega- 
tive thereby improving the longitudinal response considerably as 
seen in Fig. 2. 28. The fluctuations of M^and M^just described are 
the direct result of horizontal tail stall and subsequent blanking in 
this region of angle -of-attack as discussed in Section 2. 1. 

The longitudinal motion of the aircraft is slowly divergent 
at very large angles -of-attack for the practical range of airspeeds. 

The trend of stability with airspeed variation is reversed once again 
in that the instability becomes more severe as airspeed is decreased. 

In general, it is concluded that at low and very large angles - 
of-attack for all practical airspeeds of the transition flight regime the 
rate of divergence is slow enough that a pilot could maintain longi- 
tudinal control of the aircraft. However, in the range l6°^0< w <38 0 
the rate of divergence is such that the pilot could maintain control 
only at a reduced airspeed. If a completely flexible airspeed angle- 
of-attack capability is required or desired an automatic stabilization 
system would be necessary in the longitudinal control loop. 

Because of the limitation of time, an analysis of the sta- 
bility in the negative angle -of-attack range was not carried out. 
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However, some qualitative conclusions can be drawn about stability 
in this range from the observations of stability in the positive angle - 
of-attack range. At low and very large negative angles -of-attack 
the longitudinal stability characteristics of the aircraft could be 
expected to be similar to those already discussed for positive angles - 
of-attack. The stall would occur at a smaller negative angle than 
that positive angle at which it occurs but this would not alter the 
nature of the dynamic response of the aircraft appreciably. The 
rapid deterioration of stability in the deep stall region of positive 
angle -of-attack would not occur in the negative range since the hori- 
zontal tail would never be blanked by separated flow from the wing. 

For a number of reference flight conditions for which it 
was felt that a significant effect might be noticed, the nose fan pitch 
damping was removed from the analog circuit. Fig. 2.28 shows 
that a very slight reduction in damping resulted. There was notice- 
able change in period of oscillation with the removal of nose fan 
damping. 

Main fan lift damping, although included in the analog 
circuit, was known to be small in comparison to the wing damping 
and would not change the basic characteristics of the aircraft. 
Therefore computer runs were not performed with fan lift damping 
removed. 

2. 3.2 Lateral-Directional Mode 
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The periods and frequencies of oscillations of the lateral- 
directional responses are presented in Table II as determined from 
the analog recordings of Figs. 5. 35 through 5. 57. The period of the 
lateral-directional oscillatory response varies directly with angle- 
of-attack and inversely with airspeed. This dutch roll motion ex- 
hibited periods ranging from 3. 75 seconds at V Q = 100 knots and 
= 8° to 20 seconds at V n = 20 knots at = 72°. Considering 
the periods of the longitudinal responses and those of the lateral- 
directional mode above, it can be seen that all periods are consider- 
ably larger than one second. It is concluded, then, that pilot induced 
oscillations caused by coupling of control deflections with the natural 
frequency of the oscillations would not constitute a problem in the 
transition flight regime for this aircraft. 

The lateral-directional damping characteristics of the air- 
craft over the ranges of speed and angle -of-attack are shown in Fig. 
2.29 by means of the l/T parameter as described in the previous 
section. In this case it is clear that there is a deterioration of sta- 
bility with reduced speed in most of the angle -of-attack range. Ex- 
cept for a small range of angles -of-attack just above the stall, the 
degree of this deterioration is not as severe as it is in the longitud- 
inal mode. It is noted that, generally, the aircraft will not become 
unstable in the lateral-directional mode until speed is reduced to the 
vicinity of 40 knots. 

Within the small range of angle -of-attack, 16° to 19°, the 
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lateral-directional motion becomes rapidly divergent because of the 



Comparison of Figs. 2. 28 and 2. 29 indicates that the peak rate of 

lateral-directional divergence for V Q - 60 knots is approximately 

double that of the longitudinal divergence at the same speed. 

An interesting conclusion drawn from Fig. 2.29 is that 

blanking of the vertical tail in the separated flow behind the stalled 

wing, that is in the angle-of-attack range 24°^o< <40°, seems to 

w 

have a serious destabilizing effect only at higher speeds. At 20 knots 
blanking of the vertical tail appears to have little or no effect on the 
lateral-directional response. 



cerning the negative angle-of-attack region as far as the lateral-dir- 
ectional motion is concerned. Blanking of the vertical tail would not 
occur and hence no deterioration of stability would be experienced due 
to this effect. However, the rapid divergence in a narrow angle-of- 
attack band immediately following the stall in the negative angle region 
could still be expected to occur. 



log circuit for some runs for which a significant effect might be ex- 
pected. In only one case, at <^= 28° and V Q = 60 knots, did the fan 
damping prove to be important. Under these conditions the time to 
damp to half amplitude of the dutch roll oscillation was increased 
by a factor of three. 



change in polarity of the impicjcrtant 




Conclusions similar to the longitudinal mode are drawn con- 



Main fan roll and yaw damping were removed from the ana- 
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It is evident from the analog traces of Figs. 5. 35 through 
5. 57 that the aircraft exhibits, a very slow spiral divergence. This 
minor instability was not rapid enough to be of interest. 

From the observed lateral-directional stability charac- 
teristics it is concluded that for flight -at any transition speed over 
the positive and negative angles -of-attack ranges, except zfcl6°fr 
cK 2 0° , the rates of divergence are such that a pilot could 
maintain lateral-directional control with a manual control sytem. 
Within the range il6°^ O< w <i20° the very rapid rate of divergence 
of the dutch roll oscillation would prevent a pilot from maintaining 
control with a manual control system. Therefore, sustained flight 
in this range of angle -of-attack would make necessary automatic 
stabilization of the lateral-directional Control loops. It. is suggested, 
however, that a pilot could transit through this critical range with 
manual controls provided he remained within this angle -of-attack 
range for only a short period of time. 



The dutch roll parameter , J J, which indicates the rela- 
tive magnitude of the roll motion to the yaw motion, was measured 
directly from the analog traces. The results are plotted as a func- 
tion of angle -of-attack in Fig. 2. 30. At low angles -of-attack | 
increases as speed decreases due to the deterioration of roll damp- 
ing with reduced airspeed. For angles just above the stall J 
becomes much greater than unity. Noting, from Fig. 2.29, that 
the aircraft is very unstable in this range of angle -of-attack, it is 
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concluded that the instability is primarily a roll divergence. For 



speed. At 60 knots in this angle range, where the aircraft is stable, 
the rolling motion will be equal to or greater than the yawing motion 
whereas at 20 knots, where the aircraft is unstable, the yawing mo- 
tion will be more pronounced than the rolling motion. 

As a final indication of the regions of transition flight 
which would be hazardous with manual controls only, Fig. 2. 31 is 
presented. This figure, in which relative airspeed is plotted against 
angle -of-attack with time to double amplitude, as a parameter, 

is a cross plot of the curves of Figs. 2. 28 and 2. 29. Sustained 
flight within the speed-angle -of-attack regions outlined by the dashed 
boundaries is considered undesirable if manual controls alone, are used. 
The particular boundary used to define this prohibitory flight regime 
is dependent upon the minimum value of Time to Double Amplitude, 

T£> with which a pilot can operate. It can be seen that the lateral- 
directional boundary is very narrow in comparison to the longitudinal 
boundary although the rate of divergence is faster within the former, 
as previously discussed. The boundaries are approximate in nature 
and are not meant to indicate exact areas. 



large angles-of-attack, Fig. 2. 30 shows 




with air- 
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3.0 CONCLUSIONS 



1. Lift and drag forces at large angles of attack result from 'flat 
plate’ type forces with no contribution from circulation effects. 

2. Blanking of the horizontal tail causes rapid deterioration of static 
and dynamic longitudinal stability in the range + 20°^o< w < + 40°. 
Blanking of the vertical tail causes a moderate deterioration of lateral- 
directional dynamic stability in the range + 24°tSo< w ^+ 40° at the 
higher transition speeds. 

3. Longitudinal stability derivatives can be considered constant for 
any given reference angle -of-attack. 

4. Periods of longitudinal and lateral-directional oscillations over 
the ranges of speed and angle-of-attack considered are, in all cases, 
greater than 2. 5 seconds. Therefore, the danger of pilot induced 
oscillations in the transition flight regime is not considered to be a 
problem for the XV- 5A. 

5. At low and very large angles-of-attack there is a moderate decrease 
in longitudinal stability with airspeed reduction. In the range + 16° < 
°< w — + 36° the aircraft exhibits a severe longitudinal instability, the 
rate of divergence varying directly with airspeed. The shortest ob- 
served time to double amplitude was 0. 5 seconds at o< =? +28° and V 0 = 

w 

100 knots. 

6. At low and very large angles-of-attack there is a moderate decrease 
in lateral-directional stability with reduction of airspeed. In the range 
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+ 169=5 cx w ^ +24° the aircraft exhibits a very rapid roll divergence. 
The shortest observed time to double amplitude was 0.43 seconds at 
= + 18° and V 0 = 60 knots. 

7. Except for the range of angle -of-attack + l6°£of w £ + 40° rates 
of divergence in regions of instability are such that a pilot could 
maintain manual control about all three axes. For sustained flight 
in the excepted range, automatic stabilization of the aircraft in all 
control loops would be required. 

8. In general, the contributions of nose fan pitch damping and main 
fan roll and yaw damping to the overall damping characteristics of 
the aircraft are small. 

9. The variation of the Dutch Roll Parameter with angle -of-attack 
indicates that the lateral-directional instability in the range + 16°< 
0< w S+ 24° is primarily one of roll divergence. The maximum 
observed value of I ^Ar e l was 2. 32 at <= s <^= + 18° and V Q = 60 knots. 



34 



4. 0 BIBLIOGRAPHY 



1. Andrews, E. J. Notes on Simplified Stability Derivatives . 
U. S. Naval Postgraduate School. Monterey, California. 
Unpublished. 

2. Douglas Aircraft Company, Inc. USAF Stability and Control 
DATCOM by D. E. Ellison and L. V. Malthan. October, 
I960, revised July 1963. Contract AF 33(6l6)-6460. 

Project 8219 

3. Etkin, B. Dynamics of Flight. John Wiley and Sons, Inc. 
1959 

4. Johnson, C. L. Analog Computer Techniques. McGraw- 
Hill Book Co. 1963 

5. McCracken, Daniel D. A Guide to FORTRAN Programming. 
John Wiley and Sons, Inc. 1949 

6. National Advisory Committee For Aeronautics. Approx- 

w '* 

imate Relations and Charts for Low-speed Stability Deriva- 
tives of Swept Wings, by Thomas A. Toll and M. J. Queijo. 
Technical Report No. 1581. 1948 

7. National Advisory Committee for Aeronautics, Langley 
Aeronautical Laboratory. Effect of an Unswept Wing on the . 
Contribution of Unswept Tail Configurations to the Low-speed 
Static and Rolling Stability Derivatives of a Midwing Airplane 
Model, by William Letko and Donald R. Riley. Technical 



3'5 



Note 2175. August, 1950 

8. National Advisory Committee for Aeronautics, Langley 
Aeronautical Laboratory. Wind Tunnel Investigation of 
Effect of Sweep on Rolling Derivatives at Angles of Attack 
up to 13° and at High Subsonic Mach Numbers, including 

a Semiempirical Method of Estimating the Rolling Deri- 
vatives, by James W. Wiggins. Research Memorandum 
L54C26. 2 7 April 1954 

9. Northrop Aircraft, Inc. Dynamics of the Airplane. 
September, 1952. Bureau of Aeronautics, Navy Depart- 
ment. Report AE-61-4II. 

10. Perkins, Courtland D. and Hage, Robert E. Airplane 
Performance, Stability and Control. John Wiley and 
Sons, Inc. 1949 

1 1 . Royal Aeronautical Society. Data sheet Aircraft 00. 00. 06 
Conversion of Stability Derivatives for a Change of Body 
Axes. October 1954. 

12. Ryan Aeronautical Company. Estimated Static Stability 
and Control Characteristics of the U.S. Army XV-5A 

Lift Fan Research Aircraft. Report No. 64B031. 6 March 

1964 

13. Ryan Aeronautical Company. Estimated Dynamic Sta- 

bility Characteristics of the U.S. Army XV-5A Lift Fan 
Research Aircraft. Report No. 64B104. 18 August 1964 



36 



TABLE I 



Periods and Frequencies 
Longitudinal Mode 



Run 


Period 

Sec- 


U c 

Rad/Sec 


■ LO 1 


39-0 


0-161 


LO 2 


390 


0-161 


LO 3 


38-0 


0-165 


L04 


200 


0-314 


L04a 


20-0 


0-314 


LO 5 ' 


— 


• 


L 0 6 


— 


— 


LO 7 


— 


— 


LO 7ci 


— 


— 


LC 8 


2-65 


2-370 


L09 


. 3 -40 


2 -GOG 


L0 10 


— 


— 


LO1 1 


— 




L012 


— 





L013 


— 


— 


L013a 


— 


— 


L0 14 


— 


— 


L015 


— 


— 


L0 1 6 


34-1 


0-184 


L0 1 7 


5 1 -2 


0-103 


L018 


t 20-1 


0-313 


L019 


. 


— 


L019a 


. 


— 



Dashes signify no oscil laticr, recorded 
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TABLE II 

Periods and Frequencies 
La t e ra I- D i rectional Mode 



Run 


Period 

Sec 


0) c 

Rad / Sec 


LD 1 


3-75 


1-675 


LDIa 


3-80 


1-650 


LD 2 


5- 80 


1-081 


LD2a 


5-70 


1-100 


LD 3 


O 

O 


5-710 


LD 4 


5-20 


1-210 


LD 4a 


5 -00 


1-250 


LD 5 


980 


•640 


LD6 


— 


— 


L D 7 


— 


— 


LD8 


8-70 


•772 


L D 8a 


8- 60 


• 733 


L D 9 


13-00 


• 483 


LDS a 


12 -50 


•500 


LD 1 0 


9-50 


• 6 60 


LDlOa 


9- 20 


• 6 80 


LD11 


16-00 


• 392 


LD 1 1 a 


1 5-90 


• 394 


LD 12 


10-50 


•599 


LD 1 3 


18-50 


•339. 


LD14 


' 7-00 


' -896 


LD15 


15-50 


• 405 


LD1 6 


20-00 


• 314 



Dashes signify no osci nation recorde d 
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5. 0 APPENDIX 



5, 1 Aircraft Equations of Motion 

The equations of motion used to define the dynamic behavior 
of the aircraft after being preturbed from a reference flight condition 
are those commonly used in a dynamic stability analysis. The devel- 
opment of these equations is presented in detail in many references 
and will not be shown here. (3, 9, 10) Only the final equations and 
the assumptions upon which they are based are presented. The sym- 
bology used is that of Etkin. (3) The analysis considers only stick 
fixed conditions; therefore no contributions arising from control 
deflections are included. 

The equations of motion were developed utilizing an Eulerian 
set of axes fixed to, and moving with, the aircraft mass center. For 
ease of definition of motions along and about these axes, a body axis 
system, as defined by Etkin, was chosen. (3) Fig. 5. 1 presents the 
axis system and defines positive directions of forces, moments, and 
motions. The assumptions made in defining the equations are: 

1. The aircraft mass is considered constant. 

2. The aircraft is considered a rigid body. No aeroelastic 
effects are included. 

3. Gyroscopic effects of spinning masses are neglected. 

4. Higher order acceleration terms are neglected. 

5. The motion of the aircraft consists of small deviations from 
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some reference flight condition. 

6. The disturbance quantities and their derivatives are small 
so that squares and products are negligible in comparison to first 
order quantities. 

7. The y plane is a plane of symmetry. Products of inertia 

involving the / Y ^ and directions are therefore identically equal to 

zero. 

The reference flight condition chosen was that of symmetric 
steady state motion. This signifies the absence of any steady state 
linear acceleration and rotational velocity or acceleration. Also 
inherent in this reference flight condition is the fact that V~ a , 4^ > 

and are identically equal to zero. 

Based on the aforementioned conditions and assumptions, the 
six equations of motion may be divided into two sets of equations 
defining motions which are independent of each other. These equa- 
tions may be expressed as follows. 

Longitudinal Mode 

A X - (Knj COS 0j0 = m ( u ■+ oSl 

A2 -(my Sin So) © = 

AM = 



(5-1) 

(5-2) 

(5-3) 



Lateral-Directional Mode 
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AY + m^C^cosGo + ^ Sin ©J = m (v- +u 4 r -uj- o p) (5-4) 



All 


= A p 


- E r 


(5-5) 


AN 


= C r 


-Ep 


(5-6) 



The force and moment terms AX, AY » A 2 » All, An, and 
AH consist of a sum of contributions from various aerodynamic 
effects only. 

Generally the above equations are completely non-dimension- 
alized permitting the previously mentioned forces and moments to be 
expressed in terms of non-dimensional stability derivatives. How- 
ever, this analysis was to be conducted with the aid of an analog 
computer using dimensional derivatives throughout the simulation. 
Hence, dimensional derivatives were defined as shown in Section 5. 3 
and the final equations of motion used are as follows. 

Longitudinal Mode 

- (#r>g cos 0 O )©= m (u.+ ur 0 c£\ (5-7) 

+ +Z^u^ rr)(l>- U 0 <g) (5-8) 

> •+ ^ 4 Ah*.* = 6^ (5-9) 

Lateral-Directional Mode 

y^tr-j-Ypp ^Y r r +■ mg( <{>cos ©„ + ()* sin© 0 N )=ho(V'4Lt v ,r- ur a f) (5-10) 
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(5-11) 



L'^V- + L! p p 4- L ! r r +^p L P P = Ap - J= r 
N^. +r- 4 Hpp 4Nl r r =Cr-Ep (5-12) 

As a final simplification, engine and fan RPM was taken as 
being invariant during a perturbation in forward or vertical speed. 
The fan damping terms arise primarily from a change in air mass 
flow through the fans as the aircraft is perturbed from its reference 
flight condition. 
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Definition of Axis System with Positive Directions Shown 



5.2 Simplified Relations for some Non-dimensional Stability 



Derivatives 



In Section 2.2 it is stated that some of the estimates of stability 
derivatives were based on simplified analyses. The derivation of 
expressions for these derivatives, obtained from an unpublished ref- 
erence, are shown in this section. (1) 



5.2.1 C 



L <U ± . 



Standard expressions for and £ may be written: 



“ " L w ~ £ + L t 



= ^ £ 



( 






w 



- 



f Ai ) 



The quantity is the time required for the flow to travel from the 
wing to the tail and may be expressed: 

At = At 

u 



Then: o<\ «= ( i — \ + c) £ A c*r It _ ^ 

+ w ^ ^ d± U 



w 



L t 



From the above expression it can be seen that Ac* due to downwash 



lag is: 



A °< = 




d<x h 









alt 0 




A L = 


Tf 


7+ 




* C 4 = 


Z 


CL «t 7* 


k 



2. ^ € it c*’ c 

^°< C zo 

Z kt i±. 

C 2.U 

It 2LQ. 

^c< C ZU 
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die 



Finally: ^"L. = — - 

z u 



2 C n it it 

^ It Sw C ^ 



°<t 



(5-13) 



5.2.2 



'm. 



This derivative results directly from the product of the left 



on the horizontal tail due to downwash and the tail arm. 






ft 



= ~ Cl ^ t 



(5-14) 



5.2. 3 



' L( bt 



The derivative Ci represents the lift on the horizontal tail 

L %t 

caused by the change in angle -of-attack at the tail due to a pitch 
rate. The development is as follows. 



A <=x A 



= 

u 



* L " Z/° S t % C L« t 4^* 

*C. = V ^ S * n C, zjt qc 

L U Sw l + Lo< t c ZU 



Finally: ^ L = - 2. ^ 









°<t 



n 4^ — 

(t S w c 



(5-15) 



5.2.4 C 



*UL 



Referring to Fig. 5.2, the following expression maybe 



written. 
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X 



*T COS A 



D 



Assuming T invariant with U and differentiating with respect to U. 
gives : 



= 

u. 



c 



Xu. 



Therefore : 



_ ixiSy ( U+ a) 2 

^ a 2 / 

- C D/ >SU 

C x = -2.C 0 



\ 

TTW 



(5-16) 



5.2.5 C 






Again referring to Fig, 5.2, the following development can 
be made. 

z = - T sirtA - L 



_ q _ ^ J_ /> ( U z + 2Uu +■ u z ) C L 

^ U £ U 2 / 



_ - c, ^ US 



Z 3 



c 



= ^ c. 



= — 






,\\ 



C L /O S (j ^ 

‘ ^ S uV 



Then: 



a. - - 2c ‘ 



(5-17) 
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5. 2. 6 



From Fig. 5. 3, the following relation may be written: 

d - D = C °w 2 Z 3 ^ U ~ c ^ if 

Then: |l! = - C O w f U C d % 

Assuming a constant spanwise value of C DuJ gives: 





= - c ^/° u 8^ d f 




8 r 


= - c °«f u L^ c f d t 




c ~ 


= ^ = ( Q 2 u\ 




n r 


'Mife) " b) 






_ \ r N 2U] 






l b J 






' k 

- - C O w /’ U J 

' J ~ t/a. 


2 U 




— — ^ ^-£> w lx* ~ 


-fe ct’ 




5W 2 - 


SU'- 



Finally: C. = ~ < O. 8 (5-18) 

r «, 3 

The factor 0. 8 improves the approximation for the area moment of 
inertia, IxX* 

5 . 2 . 7 r 

The vertical tail contribution can be expressed as: 
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U z S v 

U 1 Sw 




Fig* 5*4 defines the vertical tail moment arm. The angle of sideslip 
at the vertical tail may be expressed: 





A- 


rfv - Zh . r\o 

U b Z U 




Then: 


C *Vy 


- C* 2_fy 
Av b 


(5-19) 


5.2. 8 


% 

r w 







Referring once again to Fig. 5. 3, the following relations may 
be written: 



d L - C L '/ x ( u_ r^c ciy- 

= C L I/ Z f (- 2-Uy) cdy 

^ r 

XUih » c L/ °u fed.} 



Assuming constant spanwise lift distribution gives: 



uJ 

i r 




L'r - C L /■ U I« 





C, /O U IxK 
Vyf U l 5b 



2 U 
b 
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Finally: 



(5-20) 



Cm = t Ql La 

^rw S b* 

C n — Ci_ ^ 0.8 

T 

The factor 0. 8 is again used to improve the area moment of inertia 
estimation. 

s - 2 - 9 c '’ Pw 

Fig. 5. 5 presents a spanwise element of the right wing under- 
going a positive roll rate. From this figure, the following develop- 
ment is possible. 



d M * ^{(D+AD^CoSAc* - o- ( L+AL)sm A cxL^ 

= ||^ad-lao<| 
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A. o » 8 



(5-2(1) 



Finally: 








- O 



■w 



5.2.10 Cjf 

r w 

From Fig. 5. 5, the wing contribution to roll damping may also 
be analyzed. 



d L! = - y [( L+aO C oS A»f - L 4 (D + ^0) Sin Ao<] 



= - cj [ aL 4- D A o< ] 






( C L + C d')w 




& t 







2. 

Sb r 





= cc L ^ + c D ) w r 



•X/. 






This gives: 




(5-22) 
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FIGURE 5-2 

Forces Due to Perturbati on ,u 




Spanwise VeJocity Distribution 
Due to Yaw Rate 
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FIGURE 5-4 

Vertical Tail Lever Arms 



L*AL 




FIGURE 55 

Incremental Forces on Spanwise 
Element Due to the 
'Perturb at ion, p 
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5. 3 Dimensional Derivatives - Body Axes 



The analysis of the dynamic behavior of the aircraft after be- 
ing perturbed from a reference flight condition was accomplished 
through analog simulation using a body axis system. Dimensional 
parameters were used throughout the simulation. However, initial 
analysis of the variation of the non-dimensional stability derivatives 
with angle-of-attack was conducted and presented with reference to 
a stability axis system. The method used to convert this data to 
that used for the simulation is presented in this section. 

The non-dimensional stability derivatives based on a sta- 
bility axis system were dimensionalized in a slightly different form 
than is normally used. (9) Mass and moments of inertia quantities 
were not included in order that a better physical feel for the forces 
and moments involved could be gained. The definitions of all dim- 
ensional stability derivatives used are presented in Table III. The 
calculation of these parmeters was programmed on a Control Data 
Corporation Digital Computer, Model 1604. The FORTRAN pro- 
gram used for this purpose is shown in Fig. 5.6 and, except for minor 
variations due to local procedures, is of standard format. (5) The 
input data consists of angle-of-attack and non-dimensional stability 
derivatives, in stability axes, inserted in the form of one-dimension- 
al , 23 element arrays corresponding to angles -of-attack from 0° to 
80° in 4° increments with angles -of-attack of 14° and 18° inserted 
in sequence. 
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The FORTRAN program of Fig. 5.6 was also designed to 
transform the dimensional derivatives based on stability axes to cor- 
responding derivatives based on body axes. The transform equations 
used were those contained in the RAeS Data Sheets and are presented 
in Table IV. (11) The angle of rotation, £ , through which the axes 

were rotated was taken to be the angle-of-attack of the aircraft since 
the plane containing the MAC of the wing is within 0. 5° of coincidence 
with the plane containing the body longitudinal axis. (12) 

Tables V, VI, and VII present the dimensional stability der- 
ivatives and associated angles -of-attack in body axes for relative 
velocities of 100, 60, and 20 knots respectively. 

Several derivatives which are peculiar to the particular air- 
craft used in the investigation also require definition. These der- 
ivatives are four in number and consist of the nose fan pitch damping, 
nose and main fans lift damping, and main fans roll and yaw damping. 
These contributions were derived analytically and experimentally by 
RAC and only a brief discussion of each is presented here. (13) 

Nose Fan Pitch Damping 

From wind tunnel tests, RAC found that the nose fan pitch- 
ing moment contribution could be expressed as follows. 



M 



N.F. 



AC 



N.F. 



q? s % 

UH.F n.f N-f=. 



where : 



(5-23) 
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and: 



iC v..,= s -°v,. r „ T . 

< '~t ] - V 






Expressing o< as the sum of a reference flight angle-of-attack 

* f TbT 

and a perturbation angle-of-attack, expanding in a trigonometric 
identity, and differentiating Equation (5-2 3) 

with respect to o< w/? yields: 

(- (5 ‘ 24) 



Using a small perturbation assumption, this may be expressed: 



M. 



— L.2S 






% 






'V*.* <■ *. 






— c* 



Nf= 



S i A ) 

NIP 



(5-25) 



At low angles -of-attack, the termor Sin of is very much smaller 
than coiof since both the perturbation, o< , and Sirw 0j terms are 

O/vpr "F wp 

small. This term may then be neglected in the low angle-of-attack 
region. At high reference angles-of-attack, the perturbation, 
approaches zero. This may be seen from the expression for the 

jaw 

angle-of-attack perturbation at the nose fan due to a pitch rate on 
the aircraft as derived in Fig. 5. 7 and given as: 



- - 6 CoS o< 0 (5-26) 

v 

At high angles-of-attack, then, the termo^Smo^ approaches zero. 
Therefore, throughout the angle-of-attack range 0° ^ 90° the 

derivative ^may be approximated by: 
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= 6 * z5 % - "K 




(5-27) 



The contribution of a perturbation in vertical velocity to the 



angle -of-attack change of the nose fan has been neglected in consid- 
ering pitch damping since it is very small. 

Nose and Main Fans Lift Damping 

This contribution, as well as the others to follow, was ob- 
tained by RAC directly from wind tunnel data and estimates of the 
changes of axial flow velocities through the fans. RAC shows that 
for maximum fan thrust output, the following data is obtained. 



For the analog simulation in this analysis, a nominal total value 
is used. 

7L — — 1^-0 LW / F-t / Sec. 

Fart S 

Main Fans Roll Damping 

Roll damping provided by the main fans, or lift fans, results 
from the differential lift damping forces produced by the fans during 
a rolling maneuver multiplied by an appropriate moment arm. RAC 
estimated this effect from analysis of wind tunnel tests to be: 







=• 8*2.3 I k> /f-t /Sec Per -fan 




NF 



H. 4 6" LU J F-t /S «?c- 
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Lp = ~ H 5*0 Ft Lb / (La.d. / Se* Toia.1 

This value was used in the present analysis. 

Main Fans Yaw Damping 

Analysis of wind tunnel tests by RAC also yielded an estim- 
ated value of the main fans yaw damping contribution as: 

/\L = - SIS’ Ft Lb/ /sec \ 

M. ir 

This value was likewise used in this analysis. 
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TABLE HI 



Definitions of Dimensional 
Stability Derivatives 



Longitudinal 


Lateral — D i rectiona 1 


x u = .ousc D 


Y ous r 

T v - o 

2 Y p 


Zu = - gusc l 


C?USb 

v_ 2 % 


x w- 2 ^ C L 


N - GUSb C 
N v- 2 <-n ? 


7 _eus(-c, -c n ) 

4 w _ D 


v - OUSb ^ 

P" 4 C Y P 


M w - G US C C m 
w 2 t nru« 


. C U S b e 
L >- 4 C 'P 


-7 0- Q ^ C Q 

Z w- 4 C L^ 


. CUSb 2 ^ 

V 4 C "p 


o . C S c 2 r . 
M w 4 c m,S; 


V . GU Sb „ 
Y r- 4 c Y r 


„ CUSc r 

Z d-‘ 4 ^Lq 


/ . GUSbV 
L r- 4 C l r 


_ _ GUSc% 

M q - 4 C m q 


N = OUSb 2 
N r" 4 c n r 
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TABLE 



Transform Equations 
Stability Axes to Body Axes 



X ug = X u Cos e - (Xw * Z u >SineCose ♦ Z w Sin 6 

X w B = x w Cos 2 e ♦ ( X u - Z w ) S i n eCos e - Z u S i n*e 

2 u B = Z u Cose - ( z w - x u ) S i n eCo s e - X w S i n*e 

Zw B = z .w Cos2 e *< 2 U ♦ X w )S i n eCo s e <-X u Sin 2 e 

M w C os ^ 

Zw B r Z^Co s ? e 
Mw B = M ^Cos e 
^Qg : Z q Cos e 

M v M a 

Longi tudinal 



Vv B = 


Y v 






%= 


N. v 


Co S 6 


♦ L v Sine 


Lv b = 


Lv 


Cos e 


- N v S i n e 


Yp B = 


Y P 


Cos € 


- Y r S i n' e 


Np B = 


N p 


Co S 6 


- ( N r - L p 


Lp B = 


Lp 


_ 2 
COS 6 


" ( L r * N p 


Yr B = 


Y r 


Cos e 


* Y p S i n e 


Np B = 


N r 


Co s e e 


* ( L p * N p 


Lr B = 


Lr 


Co s% 


-(Np - 1- p 



)S i n e Cos e - L r S i n e 

_ . _ z 

)Sine Cose ♦ Nr S i n e 

)SineCose ♦ L p S i n 2 e 
(SineCose - NpSl n 2 e 
#** 



Lateral - D irect ional 
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5. 4 Analog Simulation 



Analog simulation of the equations of motion of the aircraft 



was conducted on a Donner Analog Computer, Model 3100, in real 
time. A complete list of the equipment used is presented in Table 
VIII. Fig. 5. 8 shows the computer equipment set-up for simulation 
of the longitudinal equations. 



coupled between the longitudinal and lateral-directional modes, 
simulation of the two modes was conducted separately and will be 
discussed in that manner. Techniques used in the analog circuit 
design are standard. (4) 

5. 4. 1 Longitudinal Mode 



motion, as previously presented, were rearranged and used in the 
following form. 



Since the equations of motion used in the analysis are un- 



For the simulation study, the ^longitudinal equations of 



a _ Xa. U- + Xuj-UJ- - <3 cos 0 O © - L 0 - 0 a 

M J 0 



(5-28) 



im nr\ 



(m — 4 2^ ur - ^ e, 



e o e-im u o <J (5-29) 




(5-30) 



where: 




The fan contributions have been kept as separate entities 
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in order to facilitate easy removal of these quantities from the 
circuit to determine the effect of fan damping on the system. 

Figs. 5. 9 and 5. 10 show the complete analog circuit for the 
longitudinal mode. For a given angle -of-attack and free stream 
velocity, the minor derivatives, X w , Z u , Z , and M , were 

Parts °*NP 

assumed constant and inserted in the circuit as simple potentio- 
meter settings. The major derivatives, M^, and 

were not assumed constant. Their variation for + 1° of angle - 
of-^attack from the reference angle -of-attack was programmed into 
five variable base, diode function generators. The derivatives 
Zi, and Z^ were approximated by proper scaling of M^and 
respectively since, as a good approximation: 

Zi = C , M- 

Z, , qM, 

C, , and , are constants of proportionality defined by the above 
relations. This, of course, is not an exact representation of Zj>. 

or Za but the approximation is considerably better than assuming 

0 

2tZr an d Z j to be constants as the aircraft deviates from the refer- 
ence flight condition. The inability to program all derivatives on 
function generators was dictated by the limited number of function 
generators available. 

The aircraft was perturbed from its reference flight 
condition by a step input or impulse of pitch attitude, © , by means 
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of the switch network shown in Fig. 5.9. 



Table IX designates the reference flight condition for each 
investigative run. Also indicated is nose and main fan damping 
status and the © perturbation used. Tables X, XI, and XII list the 
levels and potentiometer settings used in the various runs. 

In programming the function generators for the major der- 
ivatives the values of a given derivative at the reference angle-of- 
attack and + 1° from this angle were found from Tables V, VI, and 
VII of Section 5. 3 and a convenient maximum value chosen upon 
which the voltage level was calculated. This was done to permit 
the function generators to operate near, but not at, the maximum 
voltage output of 100 volts. As an example, suppose for an 8° 
reference angle-of-attack a parameter, K , varied as below. 



h of 2. 0 giving: 





k = I. z 


3 

00 

II 

* 


k = r.<? 


0 

<r 

ii 

* 


k = I. S' 


f ^ h 


, would be based on 


= 


loo Volfs Cfl 






The maximum voltage output from the function generator, 95 volts, 
would then occur at o< = 8°. 

Since the function generators were programmed with respect 
to a specific reference angle-of-attack only the change, or pertur- 
bation, of angle-of-attack information needed to be channeled to the 
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generator as the ^ variable. This perturbation, defined in a 
manner similar to the perturbation of nose fan angle -of-attack 
shown in Section 5. 3, is: 

c< * cos o< 0 (5-31) 

V 

5. 4. 2 Lateral-Directional Mode 

The rearranged lateral-directional equations of motion 
used in the analog simulation are as follows. 



V = Yir' vr + (j^? + urj) p 


+ (Vr - 

x nn 


u^\ r + ^cose 4> +■ Sid©,, 


(5-32) 


p + Lr* 

A A 




1 1 ,i 

4- LpT 4 LfWp \o 
A A 


(5-33) 


r = e £> + 
c r c 


+ 

c 


+ N_ r r +. r 

c c 


(5-34) 



Once again the fan contributions are kept separate for the afore- 
mentioned reasons. 

Fig. 5. 1 1 depicts the complete analog circuit for the lateral- 
directional mode. Since this mode was assumed uncoupled from 
the longitudinal mode, motion in the lateral-directional case will 
not change the aircraft angle -of-attack. For a given reference angle- 
of attack and flight velocity, then, all the lateral-directional deriva- 
tives are constant. Variation of these quantities with sideslip angle 
or roll angle is also negligible for small perturbations. Hence, all 
derivatives in the lateral-directional mode were inserted as simple 
potentiometer settings. The circuit as shown in Fig. 5. 11 is 
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applicable to the areas of flight in which the derivative signs are 
identical to those at oC = 8°. To accommodate changes of sign, 
the appropriate feedback lines had to be adjusted physically. 

The aircraft was perturbed from its reference flight con- 
dition by a step or impulse input of lateral velocity, 0 T , which, 
in effect, perturbed the aircraft in sideslip angle . 

Table XIII lists the reference flight condition for each in- 
vestigative run. In addition, the status of fan damping and the at 

» r 

perturbation used is also indicated. Tables XIV and XV present 
the levels and potentiometer settings used in the various runs. 
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TABLE VIII 



LIST OF ANALOG COMPUTER EQUIPMENT 



1. Donner Scientific Company Analog Computer - Model 
3100 - with associated equipment: 

a. Dual DC Amplifier, model 3101 (15) 

b. Reference Supply, Model 3180 

c. Power Supply, Model 3105 

d. Power Panel, Model 3115 

2. Donner Variable Base Diode Function Generators - 
Models 35, 3750, and 3751 

3. Donner Electronic Multipliers - Models 33, 3730, 
and 3735 with associated Control Unit, Model 3736 

4. Data Equipment Company Plotamatic x-y plotter - 
Model 600M 

5. Brusch Instruments Mark 200 8-channel pen 
recorder 
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TABLE IX 



Analog Simulation 
Reference Flight Conditions 
Longitudinal Mode 



Run 


Vo 

Knots 


°<o- 

Deg. 


©O 

Deg. 


Nose Fa n 
Pitch Damp 


Ma i n Fa n 
Da m p i n g 


© Pert- 
Deg. 


LO 1 


1 0 0 


8 


1 0 


1 N 


1 N 


2.0 stp 


LO 2 


1 0 0 


8 


0 










i 




LO 3 


6 0 


8 


1 0 










1 .0 




LO 4 


2 0 


8 














0.5 




L O 4 a 


2 0 


8 






OUT 








1 




L O 5 


1 0 0 


16 






1 N 








imp 


L O 6 


6 0 


1 6 
















i 


LO 7 


2 0 


16 






1 










St p 


L O 7a 


2 0 


16 






OUT 








1 


LOS 


.10 0 


20 






1 N 








imp 


L O 9 


6 0 


20 




















L 0 1 0 


2 0 


20 




















L'O 11 


1 0 0 


28 




















L O 12 


6 0 


28 




















L O 13 


2 0 


28 




















L O 13a 


2 0 


36 






OUT 












L O 14 


6 0 


36 






I N 












L O 15 


2 0 


36 




















L O 1 6 


6 0 


52 
















St p 


L O 1 7 


2 0 


52 
















imp 


L O 1 8 


6 0 


72 
















St p 


L O 19 


2 0 


72 














imp 


L O 1 9 a 


2 0 


72 


1 




OUT 


1 




1 
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TABLE X 
Analog Simulation 
Invariant Levels 
Longitudinal Mode 



<2°u. 


20 Volt/Ft/Sec 




X ur 


50 Volt/Ft/ Sec 




c£* w 


100 Volt /Deg 


5730 Volt/ Rad 


<?e 


10 Volt/ Deg 


573 Volt / Rad 


Je 


100 Volt/Deg/Sec 


5730 Volt/Rad/Sec 


Xdr 


100 Volt/Ft/Sec 2 




X™ ur 


0-35 Volt / Lb 






1 00 Volt/Deg 


5730 Volt/ Rad 
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TABLE XI 



Analog Simulation 
Variable Leve I s. 
Long i tudinal Mode 



Run 


J X 

Volts 
L b. 


I Z w 

Volts 

Lb. 


X M w 

Volts 


T z . 

ck w 
Volts 


W 

Volts 


X z „ 

Volts 


Volts 


X m ~ 

Volt s 


L b. Ft. 


Lb. 


Lb.Ft. 


L b. 


Lb.Ft. 


Lb. 


Ft./ Sec. 


Ft./Seq 


Ft./Sec 


Ft./Sec 2 


Ft./Sec. 2 


Rad/ Sec 


Rad. Sec. 


LO 1 


2.5 00 


.4000 


.833 3 


.4 6.70 


2.000 


.106 6 


.0063 


.34 80 


LO 2 


\ 




\ 










\ 








LO 3 


4.170 


.66 67 


1. 3 88 










.177 7 


.0140 






LO 4 


1 2. 


50 


2 000 


4.1 66 










.53 33 


.0313 






LO 4a 


1 






1 




1 
















LO 5 


1.000 


.1 8 1 8 


.3135 


3 7 


.86 


1.66 7 


.74 10 


.-00 7 4 


.34 70 


LO 6 


1.670 


. 3030 


5 208 














.0124 






LO 7 


5.000 


. 9090 


15 6 3 














.0371 






LO 7a 


\ 








1 




} 




\ 




1 




i 


y 




LO 8 


.5000 


.12 50 


. 55 60 


4 7 


.1 2 


1.818 


.0250 


.0033 


. 3490 


LO 9 


.8333 


. 2080 


.92 70 










.0416 


. 0055 






LO 1 0 


2.500 


6250 


2.780 


\ 








.12 50 


.0164 


1 




LO 1 1 


7.14 0 


1.000 


.2 857 






.3597 


.0123 


.3500 


LO 1 2 


1 1.9 0 


1.66 7 


.4762 










.5 990 


,0206 






LO 13 


3 5 


.70 


5.000 


1.42-8 










1.798 


.0617 




* 


LO 13a 








\ 


1 










1 










LO 14 


1 1 . 


,1 1 


.1333 


.9260 










1.987 


.0926 






LO 15 


3 3.33 


.4000 


2.778 










5.960 


.2778 






LO 16 


1 3.88 


.9804 


.5128 










112 6 


.0333 






LO 1 7 


41, 


.65 


2.94 1 


1.538 










3.3 79 


. 1 000 






LOIS 


18.52 


.9260 


9260 










3.3 26 


.1042 






LO 1 9 


5 5.55 


2.17 & 


2.778 










9.980 


3 1 25 






LO 19a 


1 


\ 




1 


\ 




1 




1 


r 


1 


1 
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TABLE XEI 



Potentiometer Coefficients and Settings 
and Amplifier Gains 



Analog Simulation — Longitudinal Mode 





Pot Number 




1 


2 


3 


4 


5 


6 


7 


Run 


Problem Coefficient 




9 Cos © e 


m 


W o 


x u 

m 


1 


1 


m u 

0 


Pot Setting/Amplif i e r Ga i n 


LO 1 


.5535 


2 


.0560 


1 


.0827 


1 


.0656 


1 


.0572 


1 


.8750 


2 


.2921 


10 


LO 2 


.5630 












i 






\ 








1 






i 






1 






LO 3 


.5535 






.0336 






.0492 






.0352 






.0342 






.5245 


\ 




.1754 


1 




LO 4 










.01 12 






.0164 






.011 9 






.0114 






.3 500 


1 


.5846 


1 


LO 4 a 




























* 












J 




LO 5 










.1400 






.1 62 5 






.0925 




• 


.0082 






.3850 


id 


.2834 


10 


LO 6 










.0838 






.0974 






.0555 






.0049 






.2310 


j 


.1701 


1 




LO 7 










.0280 






.0325 






.0185 






.0016 






.7 700 


1 


.5671 


1 


LO 7 a 










ii 






j 






1 
















1 




i 


1 


LO 8 










.2800 






.2017 






.0841 






.2440 






.5600 


10 


.2771 


10 


LO 9 










.1 680 






.121 1 






.0504 






.1466 






.3365 






.1663 


\ 


LO 10 










.0560 






.0405 






.0168 






.0488 






.1120 


1 




.5549 


1 


LO 11 










0196 






.2769 






.0629 






.0169 






.7000 


1 


.2604 


1 0 


LO 12 










.0118 






.1658 






.0378 






.0102 






.4199 






.1562 


1 


LO 13 










.0039 






.0551 






.0126 






-0034 






.1400 






.5200 


1 


LO 13a 
















I 


















\ 










LO 14 










.0126 






.2080 






.0299 






.0031 






.5251 


10 


.1431 


10 


LO 1 5 










.0042 






.0694 






.0099 






.0010 






.1750 


1 




.4771 


1 


LO 16 






• 




.0100 






.2789 






.014 3 






.0054 






.7139 


1 


.1089 


10 


LO 17 










.0034 






.0928 






.0047 






.0018 






.2380 






.3630 


1 


LO 18 










.0075 






.3364 






.0021 






.0018 






.7567 






.54 60 






LO 19 










.0025 






.1120 






.0007 






.0006 






.2519 






.1815 






LO 19a 


1 








1 






1 












1 


\ 




\ 


1 
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TABLE XE 
( con t i nued ) 





Pot- Number 




8 


9 


10 


11 


12 


13 


14 


Run 




F 


’roblem Coefficient 




mg sin e. 


Zu 


1 


1 


i 

B 


1 

B 


1 

B 




Pot- Setting /Amplifier Gain 


LOI 


• 9758 


1 


1528 


10 


•1006 


10 


•0075 


1 


•1892 


1 


• 5285 


2 


-Q080 


1 


L02 


00 








\ 


























\ 


1 


\ 






L03 


•9758 






91 70 


1 


























•6345 


1 


• 54 5 3 






L04 










-3055 






























•2110 






• 1817 






L04a 










1 






1 
























\ 






•1817 






105 










•2868 


10 


•1 009 






• 0092 






• 2271 






•8914 






• 2422 


10 


L06 










-1721 






























■ 5348 






• 1453 






L07 










-5736 


1 


















1 








•1783 






•4843 


1 


L07a 










\ 




















1 








1 






\ 


\ 


L08 










-4244 


10 


• 1003 






• 0074 






•2 082 






• 2014 


10 


• 1362 


10 


L09 










-2546 


\ 


























• 12 08 


{ 


•8166 


1 


LOIO 










•8488 


1 


























• 4128 


1 


• 2723 


\ 


LO1 1 










1008 


10 


1.000 


1 


00 






0-0 






• 5353 






• 2650 


10 


LOI 2 










6067 


1 


























• 3206 






• 1589 




LOI 3 










2023 






























• 1070 






• 5301 


1 


LOI 3 a 










\ 






























\ 






J 






LOI 4 










-5019 






























• 0713 






• 8175 






LOI 5 










• 1673 






























0238 






• 2725 






L016 










• 3666 






























■ 1983 






• 1476 


10 


LOI 7 










1222 






























• 0661 






• 4922 


1 


L018 










•3987 






























• 0634 






• 8175 






L019 










• 1399 






























•0211 






• 2725 






L019a 










J 






1 




1 








1 








1 




J 






\ 
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TABLE XU 



(conti nued ) 





Pot- Number 




8 


9 


10 


11 


12 


13 


14 


Run 




Problem Coefficient 




mg sin © 0 


Zu 


1 


1 


1 

B 


1 

B 


1 

B 




Pot- Setting /Amplifier Gain 


LOI 


• 9758 


1 


-1528 


10 


•1006 


10 


•0075 


1 


•1892 


1 


• 5285 


2 


•0080 


1 


LO'2 


0-0 








\ 


























1 


\ 


\ 






L03 


• 9758 






• 9170 


1 


























• 6345 


1 


• 54 5 3 






L04 










•3055 






























■2110 






• 181 7 






L04a 










1 






1 
























\ 






•1817 






L05 










•2868 


10 


•1009 






•0092 






• 2271 






■8914 






• 2422 


10 


L06 










• 1721 






























• 5348 






• 1453 


i 




L07 










•5736 


1 


























1783 






•4843 


1 


L07a 










\ 


\ 


























J 






J 


\ 


L08 










• 4244 


10 


■ 1003 






■ 0074 






•2082 






•2014 


10 


• 1362 


10 


L09 










•2546 


\ 


























• 12 08 


\ 


• 8166 


1 


L010 










•8488 


1 


























• 4128 


1 


• 2723 


\ 


LOII 










1008 


10 


1-000 


1 


00 






00 






• 5353 






• 2650 


10 


LOI 2 










• 6067 


1 


























• 3206 






• 1589 




L013 










2023 






























• 1070 






• 5301 


1 


LOI 3 a 










\ 






























{ 












L014 










■ 5019 






























• 0713 






• 8175 






LOI 5 










• 1673 






























0238 






•2725 






L016 










• 3666 






























1983 






. 1476 


10 


L017 










• 1222 






























• 0661 






• 4922 


1 


L018 










•3987 






























• 0634 






•8175 






L0 19 










•1399 






























•0211 






•2725 






L019a 










I 






1 




1 








1 








1 




i 






l 
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TABLE XII 
(co n t i nued) 





Pot Number 




15 


16 


17 


18 


19 


Run 


Problem Coefficients 






F- 


X COS o 


<o 


Sans 


© 

perturb- 


COS<*o 




B 


V 


V 


Pot Setting /Amplifier Gain 


LOI 


■1542 


10 


0923 


1 


0800 


1 


•2000 


1 


•0117 


1 


L02 


1 


1 


\ 


























L03 • 


•5562 


1 


•1538 














• 1000 






• 0195 






L04 


•0 618 






•4614 














•0500 






• 0585 






104 a 








\ 




























L05 


• 1497 


10 


• 0896 






•0798 














• 0114 






L06 


• 5400 


1 


•1493 






















• 0189 






L07 


■ 0560 






■4479 






















• 0569 






L07a 




































LOS 


•1463 


10 


•0876 






0803 














• 0111 






L09 


•5278 


1 


• 1460 






















0185 






LOIO 


•0586 






•4379 






















• 0555 






LO1 1 


•1375 






• 0823 






•0805 














• 0104 






L012 


• 4960 






•13 71 






















•0174 






LOI 3 


• 0550 






• 4114 






















•0522 






LOI 3 a 




































LOI 4 


• 4 544 






• 1257 






















•0159 






LOI 5 


• 0505 






• 3770 






















• 04 78 






L0 16 


• 3458 






•095 6 






















• 0121 






L 01 7 


•0384 






• 2870 






















• 0364 






LOI 8 


• 1736 






• 0479 






















• 00 61 






L019 


• 0193 






1440 






















• 0183 






L019a 








1 






1 
















1 
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TABLE -X.LLl. 



Analog Si mula tion 
Reference Flight Conditions 
L at e ra I - D i re ctional Mode 



Run 


V 




© e 


Lift 


Fan 


v Pert 


O. 








Roll 


/ Yaw 


Step / Imp 




Kts- 


Deg. 


Deg. 


Dampi ng 


Ft / Sec 


LD1 


100 


8 


1 0 


1 N 


50 Step 


LDIa 












OUT 


1 




LD2 


60 










1 N 


3-75 




LD2a 












OUT 






LD3 


20 










1 N 


0-1 




LD4 


60 


16 






\ 


2 0 




LD4a 


{ 










OUT 


i ' i 




LD5 


20 










1 N 


0-5 Imp 


LD6 


60 


18 










10 




LD 7 


20 


i 










\ * 




LD8 


60 


24 






' 




5 0 Step 


LD8a 


\ 










OUT 


\ 




LD9 


20 










1 N 


0-25 




LD9a 


\ 










OUT 






LD10 


60 


28 






1 N 


2-5 




LDlOa 


1 










OUT 


i 




LD11 


20 










1 N 


0-1 




LD11 a 


\ ' 










OUT 


I * 




LD12 


60 


36 






1 N 


2-5 




LD13 


20 


\ 










0-1 




LD14 


60 


52 










2-5 




LD1 5 


20 












0-1 




LD16 


\ ’ 


72 


1 








\ 1 
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TABLE XE? 

Invariant L e v e I s — Analog Simulation 
Lateral-Directional Mode 



x v 


20 Volt /Deg 


1146 Volt /Rad 




20 Volt' / Deg 


1146 Volt/Rad 


Jr 


10 Volt/Deg/Sec 


573 Volt/Rad/Sec 


. Ip 


10 Volt /Deg /Sec 


573 Volt/ Rad/Sec 


Iv 


20 Volt/ Ft /Sec 




b 


50 Volt /Deg 


2865 Volt/ Rad 


tp 


20 Volt /Deg/Sec 2 


1146 Volt/Rad/Sec 2 


•Xr 


20 Volt /Deg /Sec- 


1146 Volt /Rad/Sec 2 
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TABLE XV 



Potentiometer Coefficients and Settings 



and Amplifier Gains — Analog Simulation 
L a te ra I - Di re ct io n al Mode 



R un 


Pot Number 


3 


4 


5 


6 


8 


9 


13 


Problem Coefficient 


♦ W. 
m 


-*--a 

m 


• 3v_ 

m 


Lv 

A 


i 

Lr 

A 


4 

A 


N v ■ 

c 


Pot. Setting /Amplifier Gain 


LD1 


• 8220 


1 


• 5781 


10 


•1653 


1 


•3006 


10 


• 2699 


10 


• 3590 


10 


• 9387 


1 


LDIa 












































LD2 


• 4921 






• 3507 






• 0992 






•1804 






• 1751 






• 2154 






• 5631 






LD2 a 






















1 






1 






1 










LD3 


• 0163 






• 11 69 






• 0331 






• 6013 


1 


• 5 836 


1 


• 7180 


1 


• 1877 






LD4 


•9737 






• 3403 






• 1003 






•2749 


1C 


•2600 


10 


•1128 


10 


•5041 






LD4a 






















; 


\ 




\ 




\ 








LD5 


0 323 






• 1139 






•0334 






•9165 


1 


• 8666 


1 


•3760 


1 


• 1682 






LD6 


■ 1689 


10 


•3368 






• 0974 






• 2380 


10 


•4432 


10 


•4769 


10 


• 4654 






LD7 


• 0361 


1 


•1 1 24 






•0325 






•9433 


i 


•1477 


J 


•1589 






•1551 






LD8 


•1431 


10 


• 3232 






•0'6 22 






•1602 


10 


• 0818 


1 


•4640 






•0806 






L D8a 






















1 




















Ll?9 


•0489 


1 


•1096 






•0207 






• 5340 


1 


■0273 






• 1513 






• 0269 






Lb 9 a 




\ 






































LD 10 


•1654 


10 


•3124 






• 0551 






• 7553 






•6011 






■1516 






• 0059 






LD10Q 




\ 






































LD 11 


• 0551 


1 


• 1040 






• 0184 






• 2513 






.2004 






• 5053 


1 


• 0020 






LD 11a 




i 




























i 








LD12 


• 2079 


10 


•2863 






• 0467 






• 5013 






•4930 






• 1746 


10 


•0637 






LD 13 


•0696 


1 


•9542 


1 


• 0156 






• 1671 






• 1643 


1 




• 6486 


1 


•0212 






LD 14 


• 2782 


10 


•2178 


10 


• 0336 






• 7317 






■1398 


10 


•161 1 


10 


• 2 583 






LD15 


• 0928- 


1 


• 7259 


1 


•01 1 2 






•2439 






• 4660 


1 


•537 0 


1 


• 0861 






LD16 


.1120 


\ 


•3630 




•0046 




1 


• 1410 






•2122 


\ 


•5720 


i 


i 


•0273 


1 
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TABLE XY 
(cont i n u e<i ) 



Run 


Pot- Number 


1 4 


16 


Problem Coeff- 


N r 

C 


N p 

C 




LD 1 


• 1219 


10 


• 2512 


1 


LD la 












L D 2 


■ 7315 


1 


• 1508 






LD 2a 














LD 3 


-24 38 






•0503 






LD4 


• 6285 






• 1979 






LD 4 a 














LD 5 


.2095 






•06 60 






LD 6 


• 3896 


1 




.7270 






LD 7 


• 1265 






• 2423 






LD 8 


• 58 58 






• 0868 






LD 8a 














LD9 


• 1953 






•0289 






LD9a 














LD10 


.1565 






• 0773 






LDlOa 














LD 11 


•0522 






• 0258 






LDIIa 














LD 1 2 


• 1610 






• 0099 






LD13 


• 0551 






• 0033 






LD 14 


•2 334 






0565 






LD15 


• 0778 






•0188 






L D1 6 


• 0179 


I 




• 0069 
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2 
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LU 

cr 
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5. 5 Analog Simulation Strip Recorder Traces 



The results of the analog simulation are shown in the follow- 
ing strip recorder traces. For the longitudinal mode, perturbations 
of U- ft. / sec., ur ft. / sec., o< degrees, ©degrees, and (^degrees/ sec. 
were recorded in response to a step or impulse perturbation of Q . 

The longitudinal responses are shown in Figs. 5. 12 through 5. 34 and 
are designated Run LO 1 through LO 19a. Table IX lists the reference 
flight conditions corresponding to each run. 

For the lateral-directional mode, perturbations nrft. /sec. , 
(^degrees, p degrees/sec., degrees, and f degrees / sec. The 
lateral-directional responses are shown in Figs. 5.35 through 5.57 
and are designated Run LD 1 through LD I £> . Table XIII lists the 
reference flight conditions corresponding to each run. 
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129 







130 




131 




132 




133 




134 




135 




136 




137 




138 




139 




140 




142 




143 




144 




145 




146 




148 




149 




150 




151 




152 




153 




154 




155 



\ 





156 




157 




158 




159 




160 




161 




162 
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